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Block Spin Effective Action for Polyakov Loops in 4D SU(2) LGT * 
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Wilhelm-Klemm-Str. 9, D-48149 Miinster, Germany 

Using a variant of the IMCRG method of Gupta and Cordery, we explicitly compute majority rule block spin 
effective actions for the signs of the Polyakov loops in 4D SU(2) finite temperature lattice gauge theories. To 
the best of our knowledge, this is the first attempt to compute numerically effective actions for the Polyakov 
loop degrees of freedom in 4D SU(2). The most important observations are: f. The renormalization group flow 
at the deconfinement transition can be nicely matched with the flow of the 3D Ising model, thus confirming the 
Svetitsky-Yaffe conjecture. 2. The IMCRG simulations of the FT SU(2) model have strongly reduced critical 
slowing down. 
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1. INTRODUCTION 

D+l-dimcnsional pure LGT's with finite pe- 
riodic extension N t in the "time" direction de- 
scribe glue systems at finite temperature 1/N t . 
They undergo a deconfinement transition, sig- 
nalled by a finite expectation value of the trace 
of the Polyakov loop. 

According to the 15 years old Svetitsky-Yaffe 
conjecture the D-dimensional effective statis- 
tical model for the Polyakov loops has an action 
(Hamiltonian) with short range interactions and 
a global symmetry group given by the center of 
SU(N). If both the deconfming phase transition of 
the gauge model and of the corresponding order- 
disorder phase transition of the spin system are 
continuous, the two models belong to the same 
universality class. 

The conjecture that the deconfinement transi- 
tion of 4D SU(2) LGT belongs to the 3D Ising 
universality class is supported by analytical cal- 
culations, and numerical estimates of its critical 
indices, c.f. 

We tested the Svetitsky-Yaffe conjecture by 
comparing flows of block spin effective actions 
for the Polyakov loops with flows of the 3D Ising 
model. Approach of both flows to a single trajec- 
tory ending in a common renormalization group 
fixed point demonstrates universality on a funda- 
mental level. 



2. BLOCKING POLYAKOV LOOPS 

Our procedure to generate and analyse block 
spin effective actions for Polyakov loops consists 
of the following steps. 

• Map the SU(2) configurations U living on 
an N t x lattice to Ising configurations 
cr(U) on an Ng lattice. The Ising variables 
are given by the signs of the Polyakov loops. 

• Block the Ising configurations with the ma- 
jority rule, using cubical blocks of size Lb- 

• Compute the effective coupling constants 
using IMCRG @. 

• Generate a renormalization group flow by 
increasing the block size Lb- 

• Compare the resulting flow with that com- 
puted directly in the Ising model. 

More formally, the effective action for the signs 
of the Polyakov loops is given by 

exp[-H'(j*)] = J DUP(n,U) exp[-S g (/7)],(l) 



with 



(blocks) 



1 + fj, x > sign ^ a x{U) 
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Here, S g is the standard Wilson action for SU(2). 
In case of an even block size Lb the sum of the 
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Ising spins <r x inside a block x' can be zero. In 
that case a positive (negative) fj, x < is selected with 
probability one half. 

Unfortunately, effective Hamiltonians contain 
an infinite number of coupling constants. In 
practical calculations one has to truncate to a 
finite set of interactions. We chose to include 
in the ansatz eight 2-point couplings and six 
4-point couplings. The 2-point couplings can 
be labelled by specifying the relative position 
of the interacting spins (up to obvious symme- 
tries): our couplings K\ . . . K% then correspond 
to 001,011,111,002,012,112,022,122. The 4-point 
couplings Kq . . . K14 are defined through Fig. 1 . 
The corresponding interaction terms in the effec- 
tive Hamiltonian are denoted by S' a , a = 1 ... 14. 







< 








/ 



K9 



Kio 



Kn 





y 




/ 


y 








1 










* 




1 

















Ki 



Ki 



Ki 



Fig. 1. Definition of the 4-point-couplings. 



3. IMCRG FOR POLYAKOV LOOPS 

Translating the rationale of IMCRG || to the 
present context means that one does not simu- 
late standard SU(2) model but instead the system 
with partition function 

Y.J DU p ^ u ) ex P [~ s ^ u ) + > ( 3 ) 



where B(fi) = ^2 a K a S' a (fi) is a guess for H'(/i). 
Note the plus sign in front of H . As before, fj, 
denote the block spins, defined on blocks of size 
Lb obtained from blocking the Polyakov loops. 
The crucial observation is now that if H = H', 
i.e., if the guess of the Hamiltonian is right, then 



the /U x / decouple completely. This means in par- 
ticular that the correlations (S' a ((i)) vanish. A 
non-perfect guess can be improved by iteration of 



K a + n- 1 (S / M), 



(4) 



where n a are trivial multiplicity factors. A few 
iterations of this correction step usually suffice to 
obtain good precision for the effective couplings. 

This method, though it seems to be restricted 
to Ising type models, has several merits. One of 
them is the following: For small enough blocks, 
the autocorrelations in the simulations are dras- 
tically reduced. This follows from the fact that 
in case of perfect compensation the blocks com- 
pletely decouple and fluctuate independently. In 
the SU(2) IMCRG calculations discussed below 
this phenomenon was clearly observed [0 . 

4. MONTE CARLO RESULTS 

We started by computing the flow of the criti- 
cal 3D Ising model, using two different models. 
The standard Ising model with nearest neigh- 
bour couplings becomes critical at f3 = 0.2216544. 
A version "I^" which includes also third (cube- 
diagonal) neighbour couplings has a critical point 
at pi = 0.128003 and /3 3 = 0.051201 @. 

We then turned to FT SU(2) with N t = 2. 
We computed the effective actions from simula- 
tions on lattices consisting of 6 3 blocks of size 
Lb < 6. The simulations were performed at 
f3 g = 1.880,1.877,1.874,1.871. Best matching 
with the Ising data was achieved at j3 g = 1.877. 
The flow of Ki and K 2 for this gauge coupling 
is shown in Fig. 2, together with the Ising flows. 
The gauge data are displayed with squares, the 
Ising results are shown with bars, diamonds and 
dotted fit lines. The fits were done with a sim- 
ple power law g. A nice matching is observed 
also for the 12 couplings not shown here. Note 
that the block sizes of the various models have 
to be rescaled with respect to each other in or- 
der to yield matching. E.g., the Lb from the J3 
model have to be rescaled by a factor of 0.59 with 
respect to that of the standard Ising model. 

Let us finally show a comparison of the flow of 
the nearest neighbour coupling for different gauge 
couplings. This is shown in Fig. 3. The dashed 
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line, together with diamonds and crosses, gives 
the Ising flow. The other data (from top to bot- 
tom) correspond to g — 1.880 (triangles), 1.877 
(crosses), 1.774 (squares), and 1.771 (stars). 

The gauge block sizes are rescaled such that 
best matching is obtained for (3 g — 1.877. Within 
the given precision, also the /3 g = 1.874 data could 
be rescaled to match the Ising flow. This is, how- 
ever, clearly ruled out for f3 g — 1.880. 
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Fig. 3. Comparison of K\ flow for different /3 g . 

5. CONCLUSIONS 

IMCRG works well as a method to compute 
the effective action of Ising type degrees of free- 
dom, also in non-Ising models like 4D FT SU(2) 
LGT. The calculations could be done on work- 
stations. The Svetitsky-Yaffe conjecture is con- 
firmed in a very fundamental way by observing 
matching of the RG trajectories with those of the 
Ising model. We obtained results for Nt = 1 also, 
see B. Extension to Nt bigger than two is ex- 
pensive, because one needs larger blocks to come 
close enough to the fixed point. 
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Fig. 2. Matching of the couplings K\ and K^. 



